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I. INTRODUCTION 



A wave function at the point x is transformed as 
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In these lecture notes I give a brief pedagogical intro- 
duction to the methods used in nonperturbative inves- 
tigations of QCD and other gauge theories. The main 
attention is payed to Wilson loops, both on the lattice 
and in the continuum, which play a central role in mod- 
ern formulations of gauge theories and to the method of 
the 1/N expansion. 

For further studies of this subject I can recommend the 
textbook [l| which contains detail references. In this text 
I restrict myself only with references to a few classical 
papers. 



Lecture I 

Wilson loops 



Wilson loops are essentially phase factors in Abelian or 
non-Abelian gauge theories. Wilson loops are observable 
in quantum theory by the Aharonov-Bohm effect. Wil- 
son loops play a central role in the lattice formulation of 
gauge theories. QCD can be reformulated through the 
Wilson loops in a manifest gauge-invariant way. Ana- 
logues of the Wilson loops are extremely useful in solving 
various kinds of matrix models. 



II. PHASE FACTORS IN QED 
A. Definition and basic properties 

Abelian phase factor is defined by the formula 



U[T yx ] = c 
Under the gauge transformation 

Af,(z) ^> A M (z) + -d^a(z) . 
the Abelian phase factor transforms as 



U[T 



iM'J 



1±> e Mv) 



y{x) 



(1) 

(2) 
(3) 



tp(x) , 



(4) 



therefore the phase factor is transformed as the product 

ip{y)ip\x): 



; V(y) f\x 



(5) 



A wave function at the point x transforms like one at 
the point y after multiplication by the phase factor: 



U[T yx ] V {x) s ~ «<p(y)», 
and analogously 

^{y)u[T yx ] e V(x)" 



(6) 



(7) 



The phase factor plays the role of a parallel transporter 
in an electromagnetic field, and to compare phases of a 
wave function at points x and y, we should first make 
a parallel transport along some contour Y yx . The result 
is F-dependent except when (z) is a pure gauge (van- 
ishing field strength F^{z)). Certain subtleties occur 
for not simply connected spaces (the Aharonov-Bohm 
effect). 



B. Propagators in external field 

Let us consider a quantum particle in a classical elec- 
tromagnetic field. To introduce electromagnetic field, 9 M 
is to be replaced by the covariant derivative 



Vjj, = 8^- leA^x) 



(8) 



For the propagator we get 



G(x,y;A) = - 



dr e 2 '' 



X / Vz (t)e~^o dti l^+ ie Io dti"(t)A M (*(t)) _ 



(9) 
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The exponent is just the classical (Euclidean) action of a 
particle in an external electromagnetic field. The path- 
integral representation © for the propagator of a scalar 
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FIG. 1: Scheme of the Aharonov-Bohm experiment. 



particle in an external electromagnetic field is due to 
Feynman. 

We can alternatively rewrite Eq. ((9]) as 

G(x,y;A) = £ e - s ^y^Ir yx ^A.(z) ^ (10) 
r 

where we represented the (parametric invariant) integral 
over dt as the contour integral along the trajectory Y yx 
over 

dz M = dti^(i). (11) 

The transition amplitude of a quantum particle in a 
classical electromagnetic field is the sum over paths of 
the Abelian phase factor |T|). 



C. Aharonov— Bohm effect 

Transverse components of the electromagnetic field de- 
scribe photons. Longitudinal components are related to 
gauging the phase of a wave function, i.e. permit one to 
compare its values at different space-time points when an 
electron is placed in an external electromagnetic field. 

In quantum mechanics, the wave-function phase itself 
is unobservable. Only the phase differences are observ- 
able, e.g. via interference phenomena. The phase differ- 
ence depends on the value of the phase factor for a given 
path T yx along which the parallel transport is performed. 

The phase factors are observable in quantum theory, in 
contrast to classical theory. This is seen in the Aharonov- 
Bohm effect Q whose scheme is depicted in Fig. [T] 
Electrons do not pass inside the solenoid where the mag- 
netic field is concentrated. Nevertheless, a phase differ- 
ence arises between the electron beams passing through 
the two slits. The interference picture changes with the 
value of the electric current. 



The phase difference depends on (the real part of) 

iej + dz»A^z) -tef dz»A„(z) 

g 1 yx g 1 ysc 

= Qie^dz^A^z) = e iefd<T»»F„„ = ^eiK ^ 

where the closed contour T is composed from Y+ x and 
r~ y . It does not depend on the shape of r^,. and T~ x but 
depends only on H S — the magnetic flux through the 
solenoid. 



III. YANG-MILLS THEORIES 

Modern theories of fundamental interactions are gauge 
theories. The principle of local gauge invariance was in- 
troduced by H. Weyl for the electromagnetic interaction 
in analogy with general covariance in Einstein's theory of 
gravitation. An extension to non- Abelian gauge groups 
was given by Yang and Mills in 1954 0]. 

A crucial role in gauge theories is played by the phase 
factor which is associated with parallel transport in an 
external gauge field. The phase factors are observable in 
quantum theory, in contrast to classical theory. This is 
analogous to the Aharonov-Bohm effect for the electro- 
magnetic field. 

A. Gauge invariance 

The principle of local gauge invariance deals with the 
gauge transformation (g.t.) of a matter field ip, which is 
given by 

ip{x) ^ ip'(x) = Sl(x)i/>(x). (13) 

Here Q (x) £ G with G being a semisimple Lie group 
which is called the gauge group (G = SU(3) for QCD). 
Equation (fT3|) demonstrates that ip belongs to the fun- 
damental representation of G. 
The unitary gauge group is when 

n -1 (x) = n + (x), (14) 

while an extension to other Lie groups is straightforward. 
Then we have 

i>\x) 2±> ip'\x) = ^{x)Q\x). (15) 

In analogy with QCD, the gauge group G = SU(N) is 
usually associated with color and the proper index of tp 
is called the color index. 

The gauge transformation (fl3|) of the matter field ip 
can be compensated by a transformation of the non- 
Abelian gauge field which belongs to the adjoint rep- 
resentation of G: 

A^^A'^x) 

= n(x) A^x) n\x) + i (l{x) d^(x) . (16) 
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It is convenient to introduce the Hermitian matrix 

[Arf = 9^2A-(x)[tT, (17) 

a 

where g is the gauge coupling constant. 

The matrices [t a ] v are the generators of G (a 
1 TV 2 — 1 for SU(N)) which are normalized such that 

trt a t b = 5 ab , (18) 

where tr is the trace over the matrix indices i and j. 

Quite often another normalization of the generators 
with an extra factor of 1/2, tr t a P' = ^S ab , is used for 
historical reasons, in particular t a — a a /2 for the SU (2) 
group, where cr° are the Pauli matrices. This results in 
the redefinition of the coupling constant, g 2 = 2g 2 . 

Equation (|17| can be inverted to give 

A%{x) = -tr A„(x)t a . 



Substituting 



finites 
SA M (x) 



we obtain for an infinitesimal a: 

= Vf<*(x). 

Here 

is the covariant derivative in the adjoint representation 
of G, while 



(19) 

(20) 
(21) 
(22) 



Snip ~ lAnljj 



(23) 



is that in the fundamental representation. It is evident 
that 



VfB(x) = [V^\B(x)], (24) 

where B(x) is a matrix-valued function of x. 

The QCD action is given in the matrix notation as 

S[A,^] 

d 4 x {dfj, - lAp) ip + rmpip ■ 



■ tr T z 
2 11 fiv 



where 



(25) 
(26) 



is the (Hermitian) matrix of the non-Abelian field 
strength. 

This action is invariant under the local gauge transfor- 
mation since 

T^{x) ^ tt{x)T^{x)rf(x) (27) 



or 



5T^{x) s = -i[J 7 flu (x),a(x)] 



(28) 



for the infinitesimal gauge transformation. 

For the Abelian group G — U(l), the formulas recover 
those for QED. 



B. Non-Abelian phase factors (Wilson loops) 

To compare phases of wave functions at distinct points, 
a non-Abelian extension of the parallel transporter is 
needed. The proper extension of the Abelian formula |1}: 



(29) 



includes the symbol P of path-ordering. 

Although the matrices A^z) do not commute, the 
path-ordered exponential on the RHS of Eq. (I29|) is de- 
fined unambiguously. This is obvious after rewriting the 
phase factor in an equivalent form 



The path-ordered exponential in Eq. (|29j) can be un- 
derstood as 

T 

U[T yx ] = Hil + idtz^^A^zit))}. (31) 
t=o 

Using Eq. (fTT|) . Eq. ([31]) can also be written as 

Up yx ] = II [l + idz^(z)]. (32) 

zer ym 

If the contour T yx is discretized, then the non-Abelian 
phase factor is approximated by 



M 



lim I 



i=l 



l+l( Zi -J!i-l)% 



(33) 



which obviously reproduces (|32p in the limit — > z^. 

The non-Abelian phase factor (f29|) is an element of the 
gauge group G itself, while A^ belongs to the Lie algebra 
of G. 

Matrices are rearranged in inverse order under Hermi- 
tian conjugation: 



u ] [T yx ] = u[v xy 



(34) 



The notation T yx means the orientation of the contour 
from x to y, while T xy denotes the opposite orientation 
from y to x. These two result in opposite orders of mul- 
tiplication for the matrices in the path-ordered product. 

The phase factors obey the backtracking {zig-zag) con- 
dition 



U[T yx ] U[T xy 



1. 



(35) 



The gauge field in the discretized phase factor 
is chosen at the center of the ith interval in order to 
satisfy Eq. ([33| at finite discretization. 

Under the gauge transformation (fT6| [/[r^J trans- 
forms as 



u\r vx ] 



g-t. 



Sl(y)U[r yx ]rt(a 



(36) 
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This formula stems from the fact that 

[l + idzM^z)] ^ [l + idzM^z)] 

= n(2! + dz)[l + idx!M M («)]nt(«) (37) 

which can be proven by substituting Eq. (fH)|) , so that 
f^(z) and Q(z) cancel in the definition (|32p at the inter- 
mediate point z. 

A consequence of Eq. (|36|) is that VK 2 -); transported by 
the matrix [/[F^] to the point y, transforms under the 
gauge transformation as ijj(y): 

U[T yx ] j,(x) g ~ >(y)". (38) 

Therefore, {/[r^] is, indeed, a parallel transporter. 

It follows from these formulas that 4>(y) E^[r yx ] ip( x ) is 
gauge invariant: 

$(y)U[r yx ]Tp(x) ^ i,{y)U[T yx ]^{x) . (39) 

Another consequence of Eq. (j36|) is that the trace of the 
phase factor for a closed contour V is gauge invariant: 

trPe^^W tT p e if r dz"A„(z) _ (40) 

This is quite similar to the Abelian phase factor. 

The sufficient and necessary condition for the phase 
factor to be independent on a local variation of the path is 
the vanishing of T^ v . The formulas of this type are well- 
known in differential geometry where parallel transport 
around a small closed contour determines the curvature. 
J-f^v in Yang-Mills theory is the proper curvature in an 
internal color space while A.^ is the connection. 

A historical remark 

An analog of thephase factor was first introduced by 
H. Weyl in 1919 [4| in his attempt to describe gravi- 
tational and electromagnetic interactions of electron on 
equal footing. What he did is associated in modern lan- 
guage with the scale rather than the gauge transforma- 
tion, i.e. the vector-potential was not multiplied by i as 
in Eq. ([T]). This explains the term "gauge invariance" - 
gauging literally means fixing a scale. 

The factor of i was inserted by London in 1927 [|[ 
after creation of quantum mechanics and the recognition 
that the electromagnetic interaction corresponds to the 
freedom of choice of the phase of a wave function and not 
to a scale transformation. 

Lecture II 

Lattice gauge theories 

Lattice gauge theories were proposed by Wilson in 
1974 Q to explain quark confinement in QCD. Lattice 



gauge theory is a nonperturbative regularization of gauge 
theory and its nontrivial definition beyond perturbation 
theory. It provides a nonperturbative quantization of 
gauge fields by a lattice. Lattice gauge theories uses 
an analogy between quantum field theory and statistical 
mechanics and offers a possibility of applying nonper- 
turbative methods: the strong-coupling expansion or the 
numerical Monte Carlo method to QCD, which provide 
evidence for quark confinement. 

IV. LATTICE FORMULATION 
A. The lattice 

A lattice approximates continuous space by a discrete 
set of points. In Euclidean formulation the lattice is along 
all four coordinates, while the time is left continuous in 
Hamiltonian approach. 

The lattice is defined as a set of points of d-dimensional 
Euclidean space with coordinates 

= n^a, (41) 

where the components of the vector 

= (m,n 2 ,...,nd) (42) 

are integer numbers. The points (|4"Tj) are called the lattice 
sites. 

The dimensional constant a, which equals the distance 
between the neighboring sites, is called the lattice spac- 
ing. Dimensional quantities are measured in units of a, 
thereby setting a — 1. 

A 2d lattice with periodic boundary conditions is de- 
picted in Fig. O Sites with the same numbers are iden- 
tified. The spatial size of the lattice has L\ = 6 and 
L2 = 4. An analogous 4d lattice is called hypercubic. 

The next concepts are the link and plaquette of a lat- 
tice which are shown in Fig. [3J A link I = {x; /x} con- 
nects two neighboring sites x and x + apt, where fi is a 
unit vector along the /^-direction (jj, = 1, . . . , d). A pla- 
quette p = {x; [i, u} is the elementary square enclosed by 

1 2 3 4 5 6 1 
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FIG. 2: Two-dimensional lattice with periodic boundary con- 
ditions. 
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g 
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FIG. 3: Link (left) and plaquette (right). 





(a) 



(b) 



FIG. 4: Description of continuum configurations by lattices. 



four links in the directions /i and v. The set of four links 
which bound the plaquette p is denoted as dp. 

For an infinite lattice, the number of degrees of free- 
dom = oo (but enumerable). To limit the this number, 
the lattice has a finite size L\ x L2 x • • • x Ld in all direc- 
tions. Periodic boundary conditions are usually imposed 
to reduce finite-size effects. 



B. Matter and gauge fields on the lattice 

Matter field, say a quark field, is attributed to the 
lattice sites, therefore a continuous field ip[x) is approxi- 
mated by its values at the lattice sites 



<p(x) 



(43) 



The lattice field tp x is a good approximation of a contin- 
uous field if (x) when a -C than the characteristic size of 
a given configuration. A description of continuum field 
configuration by lattices is illustrated by Fig. [U Lattice 
(a) is "coarse" and can represent the given continuum 
field configuration very roughly, while lattice (b) is "fine" 
and its spacing is small enough. 

The gauge field is attributed to the links of the lattice: 



A^x) 



(44) 



It looks natural since a link is characterized by a coor- 
dinate and a direction - the same as A^ (x) . Sometimes 
the notation U x ^ is used as an alternative for U^x) to 
emphasize that it is attributed to links. 
The link variable U^x) can be viewed as 



Pe 



(45) 



where the integral is along the link As a — > 0, 

this yields 



U^x) 



(46) 



so that U,j, (x) is expressed via the exponential of the /ith 
component of the vector potential at the center of the 
link. 

Since the path-ordered integral in Eq. (|45[) depends on 
the orientation, links are oriented. The same link, which 
connects the points x and x+a/t, can be written either as 
{x; /j.} or as {x + afi; — The orientation is positive for 
fi > in the former case (i.e. the same as the direction 
of the coordinate axis) and is negative in the latter case. 

The link variable U^x) is assigned to links with posi- 
tive orientations. The [/-matrices which are assigned to 
links with negative orientations are given by 



U^(x + ap) = uUx). 



(47) 



This is a one-link analog of Eq. (f34|) . 

It is clear from the relation ([45]) between the lattice 
and continuum gauge variables how to construct lattice 
phase factors (Wilson loops) - to construct the contours 
from the links of the lattice. 



C. The Wilson action 

An important role in the lattice formulation is played 
by the phase factor for the simplest closed contour on the 
lattice: the (oriented) boundary of a plaquette. The pla- 
quette variable is composed from the link variables (I45[) 
as 

U(d P ) = Ul{x) U£(x + av) U v {x + afi) U M (x) . (48) 

The link variable transforms under the gauge transfor- 
mation as 



U M (x) ^ n(x + afi)Un(x)n\x), 



(49) 



where the matrix Q(x) is attributed to the lattice sites. 
This defines the lattice gauge transformation. 

The plaquette variable transforms under the lattice 
gauge transformation as 

U(d P ) ^ n(x)U(d P )rt{x). (50) 
Its trace over the color indices is gauge invariant: 



tr U{dp) 



g .t. 



trU{dp) 



(51) 



This invariance is used in constructing an action of a 
lattice gauge theory. The simplest Wilson action is 



S la t[U] = J2 



1- ^Retr U(dp) 



(52) 
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The summation is over all the elementary plaquettes of 
the lattice (i.e. over all x, /x, and v), regardless of their 
orientations. 

Since a reversal of the orientation results in complex 
conjugation: 



trU(dp) ^ txU\dp) = [trU(dp)}" 



(53) 



the action can be rewritten as the sum over oriented pla- 
quettes 



orient p 



l_± tr £/(9p) 



(54) 



As a — > 0, the lattice action ([52| becomes the action 
of continuum gauge theory. To show this, we note that 



U(dp) -» exp[m 2 T^(x)+0(a 3 )] 



(55) 



where T^x) is non-Abelian field strength (j26|) . 

In the Abelian theory, the expansion (|55|) is easily 
found from the Stokes theorem. The commutator of 
Afj,(x) and A u (x), which arises in the non-Abelian case, 
complements the field strength to the non-Abelian one, 
as is ensured by the gauge invariance. 

A transition to the continuum limit is performed as 



(56) 



Expanding the RHS of Eq. ([55} in a, we reproduce the 
continuum action 



D. The Haar measure 

The partition function of a pure lattice gauge theory 
is defined by 



U du ^ 



x) e 



-PS[U] 



(58) 



where the action is given by Eq. 

This is the analog of a partition function in statistical 
mechanics at an inverse temperature (3 given by 



N 

P = -a 

sr 



(59) 



This formula results from comparing Eq. (|57|) with the 
gauge-field part of the continuum action ([25]) . 

A subtle question is what is the measure dll fl (x) in 
Eq. ([B"8"|) . To preserve the gauge invariance at finite lattice 
spacing, the integration is over the Haar measure which 
is an invariant group measure. The invariance of the 



Haar measure under multiplication by an arbitrary group 
element from the left or from the right: 



dU = d(flU) = d{UQ'), 



(60) 



guarantees the gauge invariance of the partition func- 
tion ([581) . 

This invariance of the Haar measure is crucial for the 
Wilson formulation of lattice gauge theories. 

An explicit expression for the Haar measure can be 
written for the SU(2) gauge group, whose element is 
parametrized by a unit four- vector a^: 



U = a^I + ma 



(61) 



where a are the Pauli matrices. The Haar measure for 
SU(2) is 



dU = in«ia^W(oJ-l) 



(62) 



since det U = a 2 ,. 

The partition function (|58[) characterizes vacuum ef- 
fects in quantum theory. Physical quantities are given 
by the averages 



(f[uy 



Z-\P) f Y[dU,(x)e-^F[U],(63) 



where F[U] is a gauge-invariant functional of U^ix). 
These become the expectation values in the continuum 
theory as a — > and (3 = N/g 2 . 



V. WILSON LOOPS ON A LATTICE 

Lattice phase factors are associated with paths drawn 
on the lattice. 

To write down the phase factor on the lattice via the 
link variables, let us specify the (lattice) contour C by 
its initial point x and by the directions (some may be 
negative) of the links forming the contour 

C = {x;/ii,...,/i n }. (64) 

The lattice phase factor U (C) is given by 

U(C) = Ufj, n (x + afii H h a/tn-i) • • ■ 

x Uf, 2 (x + a/ti) (x) . (65) 

For the links with a negative direction, it is again conve- 
nient to use Eq. ([47)1 . 

A closed contour has fa + • • • + fr n = 0. The gauge 
invariant trace of the phase factor for a closed contour is 
called the Wilson loop. 

The average of the Wilson loop is determined to be 

W(C) = ^trtf(C) 
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(R,0) 



(R,T) 




(0,0) (0,T) 
FIG. 5: Rectangular Wislon loop. 



This average is often called the Wilson loop average. 

A very important role is played by the Wilson loop 
averages for rectangular contours. Such a contour lying 
in the (x, i)-plane is depicted in Fig. [5] This rectangular 
Wilson loop is of the size R x T. 

The Wilson loop average is related for T R to 
the interaction energy of the static (i.e. infinitely heavy) 
quarks, separated by a distance R, by the formula 



(67) 



It can be proved in the axial gauge A4 = 0, where 
Ui{x) = 1 so that only vertical segments contribute to 
U(R x T). Denoting 



we have 



p e i/ R d*i^i(»i,...,i) 



/ 1 



W(RxT) = / — tr*(0)* f (T) 



(68) 



(69) 



Inserting a sum over a complete set of intermediate 
states 



(70) 



we obtain 

W(R x T) 



Y.^M\n){n\*] i {T)) 

n 

^Il^ml^fe-", (71) 



where E n is the energy of the state \n). As T — > 00, only 
the ground state with the lowest energy survives in the 



sum over states and finally we find 



W(RxT) la ^ T e - E °^ T . 



(72) 



Since nothing in the derivation relies on the lattice, it 
holds for a rectangular loop in the continuum theory as 
well. 

VI. STRONG-COUPLING EXPANSION 

Lattice path integrals can be calculated either by per- 
turbation theory in g 2 or by an expansion in f3 oc l/<? 2 . 
This is called the strong-coupling expansion in analogy of 
the high-temperature expansion in statistical mechanics 
since (3 is the analog of inverse temperature. 

To perform the strong-coupling expansion, we expand 
the exponential of the lattice action in f3. Then the prob- 
lem is to calculate the integrals over unitary group: 

I t'-tJl-t n / ' : rr--- , (73) 

where the Haar measure is normalized 



dU = 1 



(74) 



The integral (|T3"|) is nonvanishing only if n = 
m (mod TV), i.e. if n = m + kN with integer k. 

For the simplest case of m = n = 1, the answer can 
easily be found by using the unitarity of U and the or- 
thogonality relation 



1 

N 



(75) 



The simplest Wilson loop average is that for the loop 
which coincides with the boundary of a plaquette. It is 
called the plaquette average 



W(dp) 



trU(dp) 



(76) 



To calculate the plaquette average to order f3, it is 
sufficient to retain only the terms 0(f3) in the expansion 
of the exponentials: 



W{dp) 



l[ dU^ (x) [l + £ Irc tr U(dp') 

X,fi p' 



-tvU(dp) 



X ifJ, p' 



o(p 2 ) 



(77) 
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dp 



FIG. 6: Boundaries of plaquettes with opposite orientations. 



FIG. 7: Filling of a loop with plaquettes. 



The group integration can then be performed by remem- 
bering that 



aU^x){U^x)]){Ut(y)tf = ±S xy 5^5t6* (78) 



at different links. 

Using this property of the group integral in Eq. (|77|) . 
we immediately see that the denominator is equal to 1 
(each link is encountered no more than once), while the 
only nonvanishing contribution in the numerator is from 
the plaquette p' , which coincides with p but has the oppo- 
site orientation as is depicted in Fig. [6] The boundaries 
of the plaquettes p and p' have opposite orientations dp 
and dp' , respectively. 

The final answer for the plaquette average is 



W(dp) 
W{dp) 




2N 2 

t 
4 



for SU(N) with N > 3 , 
for 577(2) . 



(79) 



The result for SU(2) differs by a factor of 1/2 because 
trU(dp) is real for SU(2) so that the orientation of the 
plaquettes can be ignored. 

According to Eq. ([78]). a nonvanishing result emerges 
only when plaquettes, arising from the expansion of the 
exponentials of Eq. in f3, completely cover a surface 
enclosed by the given loop C as is shown in Fig. [7] In this 
case each link is encountered twice (or never), once in the 
positive direction and once in the negative direction, so 
that all the group integrals are nonvanishing. 

The leading order in (3 corresponds to filling a minimal 
surface: 



W{C) = [W(dp)] 



A mi „(C) 



(80) 



where W(dp) is given by Eq. (|T9")) and A m j n (C) is the area 
(in units of a 2 ) of the minimal surface. For a rectangle 



(a) (b) 

FIG. 8: Lines of force between static quarks, 
the minimal surface is a piece of the plane: 



W(R x T) = [W{dp)\ 



RT 



(81) 



to the leading order in (3. 

Contribution of more complicated surfaces, do not ly- 
ing in the plane of the rectangle, is W(C) ~ ^ arca an d 
is, therefore, suppressed since their areas are larger than 
A ■ 



VII. AREA LAW AND CONFINEMENT 

The exponential dependence of the Wilson loop aver- 
age on the area of the minimal surface (as in Eq. (|80[) ) is 
called the area law. It is customarily assumed that if an 
area law holds for loops of large area in the pure SU (3) 
gauge theory then quarks are confined. In other words, 
there are no physical \in) or (out\ quark states. This is 
the essence of Wilson's confinement criterion. The argu- 
ment is that physical amplitudes (e.g., the polarization 
operator) do not have quark singularities when the Wil- 
son criterion is satisfied. 



A. Linear potential 

A justification for the Wilson criterion is based on the 
relationship (jUTj) between the Wilson loop average and 
the potential energy of interaction between static quarks. 
When the area law 



W(C) 



large C 



-KA min (C) 



(82) 



holds for large loops, the potential energy is a linear func- 
tion of the distance between the quarks: 



E(R) 



KR. 



(83) 



Lines of force between static quarks are depicted in 
Fig. [8] for (a) linear and (b) Coulomb interaction poten- 
tials. For the linear potential the lines of force are con- 
tracted into a tube, while they are distributed over the 
whole space for the Coulomb one. 

The coefficient K is called the string tension because 
the gluon field between quarks contracts to a tube or 
string, whose energy is proportional to its length. The 
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value of K is the energy of the string per unit length. 
This string is stretched with the distance between quarks 
and prevents them from moving apart to macroscopic 
distances. 

Equation (f8"Tj) gives 

K = i ln ir = >( 2J V) m 

for the string tension to the leading order of the strong- 
coupling expansion. The next orders of the strong- 
coupling expansion result in corrections in [3 to this for- 
mula. Confinement holds in the lattice gauge theory to 
any order of the strong-coupling expansion. 

B. Asymptotic scaling 




FIG. 9: String tension versus 1/g' 2 



Equation (|84|) establishes the relationship between lat- 
tice spacing a and the coupling g 2 . Let K equals its 
experimental value 



K = (400 MeV) 2 » 1 GeV/fm 



(85) 



which results from the slope of the Regge trajectory a' — 
l/2wK, The slope a' = 1 GeV~ 2 from the p - A 2 - g 
trajectory. 

The renormalizability prescribes that variations of a, 
which plays the role of a lattice cutoff, and of the bare 
charge g 2 should be made simultaneously in order that 
K does not change. 

Given Eq. (|84|) , this procedure calls for a — ► oo as 
g 2 — » oo. In other words, the lattice spacing is large in 
the strong-coupling limit, compared with 1 fm. Such a 
coarse lattice cannot describe the continuum limit and, 
in particular, the rotational symmetry. 

In order to pass to the continuum, the lattice spacing 
a should be decreased. Equation l|84p shows that a de- 
creases with decreasing g 2 . However, this formula ceases 
to be applicable in the intermediate region of g 2 ~ 1 and, 
therefore, a ~ 1 fm. 

To further decrease a, we further decrease g 2 . While 
no analytic formulas are available at intermediate values 
of g 2 , the expected relation between a and g 2 for small 
g 2 is predicted by the known two-loop Gell-Mann-Low 
function of QCD. 

For pure SU(3) Yang-Mills, Eq. 
small g 2 by 



84)) is replaced at 



For such values of g 2 , where asymptotic scaling holds, 
the lattice gauge theory has a continuum limit. 

The knowledge of the two asymptotic behaviors says 
nothing about the behavior of a 2 K in the intermediate 
region of g 2 ~ 1. There can be either a smooth transition 
between these two regimes or a phase transition. Numer- 
ical methods were introduced to study this problem. 



C. Relation to second-order phase transition 

Continuum limits of a lattice system are reached at 
the points of second-order phase transitions when the 
correlation length becomes infinite in lattice units. 

The correlation length is inversely proportional to 
Aqcd: 



A" 1 

QCD 



a exp 



dg 2 



B(g 2 ) 



(87) 



where B(g 2 ) is the Gell-Mann-Low function. The only 
chance for the RHS to diverge is to have a zero of B{g 2 ) 
at some fixed point g 2 — g 2 . Therefore, the bare coupling 
should approach the fixed-point value g 2 to describe the 
continuum. 

B(0) = for a non-Abelian gauge theory so that g 2 = 
is a fixed-point value of the coupling constant. Therefore, 
the continuum limit is associated with g 2 — > 0. 



K = 



1 

const • 

or 



8tt 2 
llg 2 



-877 2 /llff 2 



(86) 



where the two-loop Gell-Mann-Low function is used. 

The exponential dependence of K on 1/g 2 is called 
asymptotic scaling. Asymptotic scaling sets in for some 
value of 1/g 2 . The dependence of the string tension on 
1/g 2 is shown in Fig. [9] The strong-coupling formula 
(|84[) holds for small 1/g 2 . The asymptotic-scaling for- 
mula (|B"6"|) sets in for large 1/g 2 . Both formulas are not 
applicable in the intermediate region 1/g 2 ~ 1. 



Lecture III 

1/N Expansion 

An effective coupling constant of QCD becomes large at 
large distances, so fluctuations of scales of different orders 
of magnitude are essential and there is no small param- 
eter, 't Hooft proposed in 1974 0] to use the number of 
colors N of the gauge group SU(N) as such a parameter 
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and to perform an expansion in 1/N. The motivation 
was the 1/A-expansion in statistical mechanics. 

The 1/A-expansion of QCD rearranges perturbation 
theory in a way consistent with a string picture. The 
accuracy of large- N QCD is of the order of the ratios of 
meson widths to their masses (10-15%). While QCD is 
simplified in the large- N limit, it is not yet solved. 



VIII. INDEX OR RIBBON GRAPHS 

In order to describe the l/W-expansion of QCD, it is 
convenient to use the matrix-field representation 



The matrix ([88| is Hermitian and differs from (fT7|) by g. 
The propagator of the matrix field (x) reads 

(W^)> GaUB8 = - jfF*") D^-y), 

(89) 

where we have assumed, as usual, a gauge-fixing to de- 
fine the gluon propagator in perturbation theory. For 
instance, one has 1 



(90) 



in the Feynman gauge. 

Equation (|89[) can be derived immediately from the 
standard formula 

(A«(x) 4(y)> Gauss - 6 ab D^ v (x - y) (91) 

multiplying by the generators of the SU(N) gauge group 
according to the definition ([55)1 and using the complete- 
ness condition 

N 2 -l 



(t a y j {t a ) kl = (s a s ki - —5 ij s kl 



for SU(N) 



(92) 

Alternatively, Eq. (f89|) can be derived directly from a 
path integral over the matrix fields. 

We concentrate only on the structure of diagrams in 
the index space, i.e. the space of the indices associated 
with the SU(N) group. We shall not consider, in most 
cases, space-time structures of diagrams which are pre- 
scribed by Feynman's rules. 

Omitting at large N the second term in parentheses 
on the RHS of Eq. (|89|) , we depict the propagator by the 
double line 



(A^(x)A kl (y))^ cx S il S kj 



(93) 



1 Here and in Eq. Il95|l <5 MI , is due to Euclidean metric and is to be 
substituted by —g^v in Minkowski space. 



Each line, often termed the index line, represents the 
Kronecker delta-symbol and has an orientation which is 
indicated by arrows. This notation is obviously consis- 
tent with the space-time structure of the propagator that 
describes a propagation from x to y. 

Arrows are a result of the fact that the matrix A*j is 
Hermitian and its off-diagonal components are complex 
conjugate. 

Double lines appear generically in all models describing 
matrix fields in contrast to vector (in internal symmetry 
space) fields, whose propagators are depicted by single 
lines. 

The three-gluon vertex is depicted in the double-line 
notation as 




cx g (5 l «3^2ji^3j2 _ gun^m^n^ ( ( 94 ) 

where the subscripts 1, 2 or 3 refer to each of the three 
gluons. The relative minus sign arises from the commuta- 
tor in the cubic- in- A term in the action (|25|) . The color 
part is antisymmetric under an interchange of gluons. 
The (momentum space) space-time part 



7miM2M3(pi>P2,P3) =<W 2 (P1-P2), 



P3) 



(P3-Pl), 



(95) 



is also antisymmetric. We consider all three gluons as 
incoming so their momenta obey P1+P2+P3 = 0. The full 
vertex is symmetric as is prescribed by Bose statistics. 

The four-gluon vertex involves six terms - each of them 
is depicted by a cross - which differ by interchanging of 
the color indices. We depict the color structure of the 
four-gluon vertex for simplicity in the case when i\ = 
h =i,i2 = h = h *3 = 34 =k,i 4 = ji = I, but k, I 
take on different values. Then only the following term is 
left: 



oc g* 



(96) 



where there are no delta-symbols since the color structure 
is fixed. We pick up only one color structure by equating 
indices pairwise. 

Diagrams of perturbation theory can now be com- 
pletely rewritten in the double-line notation. The sim- 
plest one describing the one-loop correction to the gluon 
propagator is depicted in Fig. [TDJ The sum over the N 
indices is associated with the closed index line. The con- 
tribution of this diagram is ~ g 2 N ~ 1. 

In order for the large- N limit to be nontrivial, the bare 
coupling constant g 2 should satisfy 



1 

N 



(97) 
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FIG. 10: Double-line representation of a one-loop diagram for 
gluon propagator. 




FIG. 11: Double-line representation of a four-loop diagram. 



This dependence on N is also prescribed by the 
asymptotic-freedom formula 



12^ 



117V In (A/Aqcd) 



(98) 



of the pure SU (N) gauge theory. 

Thus, the contribution of the diagram in Fig. [TU] is of 
order ~ g 2 N ~ 1 in the large-iV limit. 

The double lines can be viewed as bounding a piece of 
a plane. These lines represent a two-dimensional object. 
In mathematics these double-line graphs are often called 
ribbon graphs or fatgraphs. They are connected with 
Ricmann surfaces. 



Remark on the U(N) gauge group 

The double- line representation of the diagrams holds, 
strictly speaking, only for the U (N) gauge group, whose 
generators 



T A = (t a , I/Vn\ , trT A T B = 
obey the completeness condition 

N 2 



A = l, 



,N 2 



A=l 



§ il 6 kj 



for U(N) 



(99) 



(100) 



Elements of both SU (N) and U (N) can be represented 
in the form 



U 



, iB 



(101) 



where B is a general Hermitian matrix for U(N) and a 
traceless Hermitian matrix for SU (N) . 

The large- TV limit of both the U (N) and SU(N) groups 
is the same. 



IX. PLANAR AND NONPLANAR GRAPHS 

The double-line representation of perturbation-theory 
diagrams is very convenient to estimate their orders in 
1/N. Each three- or four-gluon vertex contributes a fac- 
tor of j or j 2 , respectively. Each closed index line con- 
tributes a factor of N, while g 2 ~ 1/N. 




FIG. 12: Double-line representation of a nonplanar diagram. 



A. 't Hooft topological expansion 

Let us consider a typical diagram for the gluon propa- 
gator depicted in Fig. [TTJ The sum over the N indices is 
associated with each of the four closed index lines, whose 
number is equal to the number of loops. The contribution 
of this diagram is ~ g s N 4 ~ 1. 

Diagrams of this type, which can be drawn on a sheet 
of paper without crossing any lines, are called planar di- 
agrams. For such diagrams, the addition of a loop in- 
evitably results in the addition of two three-gluon (or one 
four-gluon) vertices. A planar diagram with ni loops has 
rii closed index lines. It is of order 



rt2-loop planar diagram 



(g 2 N) 



1, (102) 



so that all planar diagrams survive in the large- N limit. 

Let us now consider a nonplanar diagram of the type 
depicted in Fig. [T^J The diagram has six three-gluon 
vertices but only one closed index line (although it has 
three loops!). The order of this diagram is ~ g 6 N ~ 
1/N 2 . 

This nonplanar diagram can be drawn without line- 
crossing on a surface with one handle (or hole) which in 
mathematics is called a torus or a surface of genus one. A 
plane is then equivalent to a sphere and has genus zero. 
A general Riemann surface with h holes has genus h. 

The above evaluations of the order of the diagrams can 
be described by the single formula 



S('i)ns-/i diagram ~- ( 



(103) 



The expansion in 1/N rearranges perturbation-theory di- 
agrams according to their topology as demonstrated in 
1974 by 't Hooft 0. It is referred to as the topological 
expansion or the genus expansion. 

Only planar diagrams associated with genus zero sur- 
vive in the large-A^ limit. The problem of summing 
the planar graphs is complicated but simpler than that 
of summing all the graphs, since the number of planar 
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graphs with no vertices grows geometrically at large uq: 

# p (rio) = n ° °f planar graphs ~ const™ , (104) 

as shown by Tuttle [|| and Koplik, Neveu, Nussinov 
while the total number of graphs grows factorially with 
no- There is no dependence in Eq. (|104|) on the number 
of external lines of a planar graph which is assumed to 
be much less than n . 

There is a big difference between the planar diagrams 
and the ladder diagrams which describe e + e~ elastic scat- 
tering in QED. For the ladder with n rungs, there are n! 
ladder diagrams, but only one of them is planar. This 
shows why the number of planar graphs is much smaller 
than the total number of graphs, most of which are non- 



planar. 



B. Topological expansion (continued) 

Equation (|103p holds, strictly speaking, only for the 
gluon propagator, while the contribution of all planar di- 
agrams to a connected n-point Green function is ~ g n2 , 
which is its natural order in 1/N. The three-gluon Green 
function is ~ g, the four-gluon one is ~ g 2 and so on. The 
contributions of all planar diagrams are of the same order 
~ 1 in the large- N limit, independently of the number of 
external lines, for the Wilson loop average 



N 



trPc 



lg § T da 



' n=0 



dx^ 



dxf 



dx u " G 



(n) 



(105) 



with 

G$..-» n { x u ...,»„) = — (tr [A^ (xi)--- A^{x n )\ ) . 

(106) 

The factor of 1/iV, which normalizes the trace, provides 
the natural normalization. The ordering along a closed 
path implies cyclic-ordering in the index space as de- 
picted in Fig. 1131 where we omit the arrows for simplicity. 
This diagram has no = 10 vertices, n± = 12 gluon prop- 
agators, ri2 = 4 closed index lines, and B = 1 boundary. 
The color indices of the external lines are contracted by 
the Kronecker delta-symbols (represented by the single 
lines) in a cyclic order. The extra factor of 1/N arises 
from the normalization. The order in 1/N of the diagram 
in Fig.rjjis ~ 1/N 2 in accord with Eq. flUg]) . 

Analogously, the color indices in Eq. (|106[) are con- 
tracted in the cyclic order. The delta-symbols, which 
contract the color indices, are depicted by the single lines. 
They can be viewed as a boundary of the diagram. The 
actual size of the boundary is not essential - it can be 




FIG. 13: Generic double-line index diagram. 



shrunk to a point. Then a bounded piece of a plane will 
be topologically equivalent to a sphere with a puncture. 
We draw planar diagrams in a plane with an extended 
boundary (boundaries) rather than in a sphere with a 
puncture (punctures). The closed boundary is associated 
with the trace over the color indices of the multi-point 
Green function 

The boundary represents the Wilson loop = a trajec- 
tory of a heavy quark in the fundamental representation. 



C. Topological expansion and quark loops 

It is easy to incorporate quarks in the topological ex- 
pansion. A quark field belongs to the fundamental repre- 
sentation of the gauge group SU(N) and its propagator 
is represented by a single line 

(tpiiij) cx 8ij = i — » — j. (107) 

The arrow indicates, as usual, the direction of propa- 
gation of a (complex) field ip. These arrows are often 
omitted for simplicity. 

The diagram for the gluon propagator which involves 
one quark loop is depicted in Fig. 114b . It involves one 
quark loop and has no closed index lines so that its order 
is ~ g 2 ~ 1/N. The diagram in Fig. [Mp is ~ g 6 N 2 ~ 
1/N analogously. 

It is evident from this consideration that quark loops 
are not accompanied by closed index lines. One should 
add a closed index line for each quark loop in order for 
a given diagram with L quark loops to have the same 
double-line representation as for pure gluon diagrams. 
Therefore, given Eq. (|103|) . diagrams with L quark loops 
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a 



(b) 



FIG. 14: Diagrams for gluon propagator which involve quark 
loop. 



associated with a given graph, is open and has B + L 
boundaries. This surface can be closed by attaching a 
cap to each boundary. The single lines then become dou- 
ble lines together with the lines of the boundary of each 
cap. We have already considered this procedure when 
deducing Eq. (fTU5|) from Eq. (fTUgj) . 

The number of faces for a closed Riemann surface con- 
structed in such a manner is n 2 + L + B, while the num- 
ber of edges and vertices are n\ and no, respectively. 
Euler's theorem states that 



are suppressed at large N by 
L quark loops ~ 



L+2-gcnus 



(108) 



The single-line representation of the quark loops is sim- 
ilar to that of the Wilson loop. Such a diagram emerges 
in gluon corrections to the vacuum expectation value of 
the quark operator 



O = 



(109) 



where the factor of 1/N is introduced to make it 0(1) in 
the large- N limit. Therefore, the external boundary can 
be viewed as a single line associated with valence quarks. 



X = 2-2/i = n 2 +L + B-n 1 +n . (113) 
Therefore the RHS of Eq. () 1 1 2[) can be rewritten as 

^y"2— m+no-B _ jy2-1h-L-2B (114) 

We have thus proven that the order in 1 /N of a generic 
graph does not depend on its order in the coupling con- 
stant and is completely expressed via the genus h and the 
number of virtual quark loops L and external boundaries 
Bby 



generic graph 



2/i+L+2(B-l) 



(115) 



For B = 1, we recover Eqs. (fTU3"]l and (fTU5)) . 



D. The proof of topological expansion 



To prove Eqs. (|103p and its quark counterpart (|108|) . 
let us consider a generic diagram in the index space which 
has n three-point vertices (either three-gluon or quark- 

gluon ones), four-gluon vertices, ni propagators (ei- 
ther gluon or quark ones), n 2 closed index lines, L virtual 
quark loops and B external boundaries. Its order in 1/N 
is 



I (3) , o (4) 

a n o + 2n o /V™ 2 

N By 



(3)/o (4) D 

7V" 2 ~"o /2-«o - B (HQ) 



as has already been explained. The extra factor of 1/N B 
arises from the extra normalization factor of 1/N in op- 
erators associated with external boundaries. 

The number of propagators and vertices are related by 



2ni 



3n, 



(3) 



4n, 



(4) 

' 



(111) 



since three- and four-point vertices emit three or four 
propagators, respectively, and each propagator connects 
two vertices. We then rewrite the RHS of (|110[) as 



N % 



(112) 



where no — + is the total number of vertices. 

The exponent on the RHS of Eq. (|112[) can be ex- 
pressed via the Euler characteristic x °f a given graph 
of genus h. An appropriate Riemann surface, which is 



E. 



't Hooft versus Veneziano limits 



In QCD there are several species or flavors of quarks 
(it-, d-, s- and so on). We denote the number of flavors 
by Nf and associate a Greek letter a or j3 with a flavor 
index of the quark field. 

The quark propagator then has the Kronecker delta- 
symbol with respect to the flavor indices in addition to 
Eq. ((Toil): 



cx S al3 S, t 



Their contraction results in 



N f 



N f . 



(116) 



(117) 



Therefore, an extra factor of Nf corresponds to each 
closed quark loop for the Nf flavors. 

The limit when Nf is fixed as N — > oo is called the 
't Hooft limit. Only valence quarks are then left (the 
quenched approximation) . In order for a meson to decay 
into other mesons built out of quarks, a quark-antiquark 
pair must be produced out of the vacuum. Consequently, 
the ratios of meson widths to their masses are 



Ttotal 

M 



N 



(118) 



in the 't Hooft limit. The ratio on the LHS of Eq. (fTTSj) 
is 10-15% experimentally for the p-meson. The hope of 
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solving QCD in the 't Hooft limit is the hope to describe 
QCD with this accuracy. 

An alternative large- N limit of QCD, when Nf ~ N 
as N — > oo, was proposed by Veneziano in 1976 [lOj. A 
general diagram with L quark loops will contribute 

/ N \ L ( 1 \ gcnus 
L quark loops ~ j ^— J , (119) 

since each quark loop results in AT/. 

The quark loops are not suppressed at large N in the 
Veneziano limit 

N f ~ N -► oo (120) 

if the diagram is planar. 

It is the Veneziano limit that is related to the hadronic 
topological expansion in the dual-resonance models. In 
the Veneziano limit hadrons can have finite widths ac- 
cording to Eq. (fTTS)) . 

X. LARGE- TV FACTORIZATION 

The vacuum expectation values of several colorless or 
white operators, which are singlets with respect to the 
gauge group, factorize in the large-iV limit of QCD (or 
other matrix models) as was first noticed by A. A. Migdal 
and independently E. Witten in late 1970's. 

The simplest gauge-invariant operators in a pure 
SU(N) gauge theory are the closed Wilson loops 

$(C) = — trPe'^c^W . (121) 

They obey the factorization property 

( • • • $(C„) ) = ( $(d) ) • • • < $(C„) ) + 0(N- 2 ) . 

(122) 

The factorization implies a semiclassical nature of the 
large- iV limit of QCD (a saddle point in the path integral 
for certain variables). 

The factorization property also holds for gauge- 
invariant operators constructed from quarks as in 
Eq. (|109p . For the case of several flavors Nf, we nor- 
malize these quark operators by 

Or = -^I>. (123) 



Here T denotes one of the combination of the 7-matrices: 
r = I, 75, 7m, nV75, = ^[7m>7^ >••• ■ (!24) 

The factorization of the gauge-invariant quark opera- 
tors holds both in the 't Hooft and Veneziano limits: 

( Or, • • • Tn ) = ( Tl ) • • • ( Vn > + 0(l/(N f N)) . 

(125) 

The nonfactorized part, which is associated with con- 
nected diagrams, is ~ 1/N in the 't Hooft limit. This 
leads, in particular, to the coupling constant of meson- 
meson interaction of order 1/N. The Veneziano limit is 
analogous to pure Yang-Mills. 

The factorization can be seen (at all orders of per- 
turbation theory) from Eq. (|1 15[) for the contribution of 
a generic connected graph of genus h with B external 
boundaries which are precisely associated with the Wil- 
son loops (or the quark operators Or)- The diagrams 
with gluon lines emitted and absorbed by the same oper- 
ator are products of diagrams having only one boundary. 
Their contribution is of order one. The diagrams with 
gluon lines emitted and absorbed by two different oper- 
ators have two boundaries. This proves the factorization 
property (|125|) at all orders of perturbation theory. 

The large- N factorization can also be verified beyond 
perturbation theory at all orders of the strong-coupling 
expansion in the SU(N) lattice gauge theory. A non- 
perturbative proof of the factorization was given using 
quantum equations of motion (the loop equations) [111 ]. 



XI. CONCLUSION 

We have considered in these lecture notes the basic 
features of the methods for nonperturbative studies of 
gauge theories, which were developed in the second half 
of 1970's - early 1980's. Their contemporary applica- 
tions in high-energy physics are extremely broad: from 
the scattering of particles at very high energies to the at- 
tempts of constructing a unified theory of all interactions, 
including gravity. Some of these issues are considered in 
other lectures at this School. 
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